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ABSTRACT 

Context. Only two planetary systems around old ms-pulsars are currently known. Young radio pulsars and radio-quiet neutron stars 
caimot he covered by the usually-applied radio pulse timing technique. However, finding suhstellar companions around these neutron 
stars would be of great interest - not only because of the companion's possible exotic formation but also due to the potential access 
to neutron star physics. 

Aims. We investigate the closest young neutron stars to search for substellar companions around them. 

Meti:)Ods. Young, thus warm substellar companions are visible in the Near Infrared while the neutron star itself is much fainter. Four 
young neutron stars are moving fast enough to enable a common proper motion search for substellar companions within few years. 
Results. For Geminga, RX J0720.4-3125, RX J1856.6-3754, and PSR J1932+1059 we did not find any co-moving companion down 
to 12, 15, 11, 42 Jupiter masses for assumed ages of 1, 1, 1, 3.1 Myrs and distances of 250, 361, 167, 361 pc, respectively. Near 
Infrared limits are presented for these four as well as five other neutron stars for which we currently have only observations at one 
epoch. 

Conclusions. We conclude that young isolated neutron stars rarely have brown dwarf companions. 

Key words, giant planet formation - neutron stars - pulsars - RX J0720.4-3125, Geminga, RX JI856.6-3754, PSR JI932-I-I059 



1. Introduction 

For a long time substellar companions (planets and brown 
dwarfs) around neutron stars have not been expected due to 
the powerful supernova giving birth to the compact object. 
Astonishingly however, the first planetary mass companions 
ever discovered around other stars were neutron star (NS) 
companions. Three such objects are confirmed around the old 
minisecond(ms)-pulsar PSR 1257-1-12 dWolszczan & Fraill 19921: 
iKonacki & Wolszczanl2003h . They have up to four Earth masses 
and are very close to the primary - 0.2 AU to around 0.5 AU. 
Another planetary system around a NS is the binary system 
of the ms-pulsar PSR B 1620-26 and its white dwarf co m- 
panion jSigurdsson et alj|2003l;lsigurdsson & Thorsettll2005h . A 
2.5 Jupiter mass object is orbiting the white dwarf at 23 AU 
and thus also the NS. Currently, we do not know of any brown 
dwarf companion around a neutron star. Several formation mod- 
els have been discussed to explain the unexpected discovery of 
NS planets either by assuming supernova survival, by planet for- 
mation after the supernova or by planet capturing (for a detailed 
review see Podsiadlowski 1995). The planet formation from in- 
situ disks around NSs could differ from those around solar-like 
stars. Thus, the study of NS planetary systems in comparison to 
those around solar-like stars can be very fruitful for the under- 
standing of planet formation in general. But substellar compan- 
ions of NSs are not only interesting for planet or brown dwarf 
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formation theory. They could also provide insights into the su- 
pernova event and maybe the composition and final stages of the 
progenitor star(s). Studying low-mass companions around NSs, 
especially the companion masses and their orbital characteris- 
tics, can determine also the dynamical masses of the NSs. These 
masses, together with the star radii, constrain the equation of 
state of matter in the neutron stars. Currently, no radius and mass 
of a neutron star are known precisely for the same object and the 
equation of state of ultra-dense matter in the universe is still not 
known. Whereas masses of N Ss can be derived for binary sys- 
tems (e.g. iJonker et al.ir2003b . the best radius estimation comes 
from the radio-quiet X-ray t hermal NSs, which ar e isolated with- 
out stellar companions (e.g.l Trtimper et all2004l) . 
Only the two planetary systems around ms-pulsars mentioned 
above were found so far, both by the radio pulse timing tech- 
nique, which measures variations in the pulse arrival times. This 
technique is most sensitive to very stable ms-pulsars and requires 
much observing time. For ordinary radio pulsars the timing noise 
hinders the precise measurements needed. Radio-quiet neutron 
stars can not be searched for substellar companions by this 
method. The radio pulse timing technique can detect only rel- 
atively close companions like in radial velocity planet searches 
among normal solar-like stars. Therefore, substellar compan- 
ions in wide orbits cannot be detected with this method. As 
NSs are very faint at visible wavelengths, the commonly used 
radial velocity method or the search for astrometric wobbling 
should preferably be done in X-rays where many NSs are bright. 
Unfortunately the spectral and spatial resolution in X-ray astron- 
omy is currently not good enough for these methods. Transit 
searches would not be effective, because rarely more than one 
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NS would be in one field of view. 

Substellar companions of NSs at any separations can be found 
only by their common proper motions detected by direct imag- 
ing, primarily in the near infrared (NIR) if the companions are 
young. They are several magnitudes fainter in the optical. The 
companions are brighter than the NSs in the NIR, even though 
lower in mass, because they are larger in size and still contract- 
ing and, hence, warm when young. Direct detection of very close 
(even few AU) and faint substellar companions is currently pos- 
sible only around NSs. Substellar companions in wide orbits can 
be detected too. However, one is restricted to young objects at 
small distances from the sun. The magnitudes of substellar ob- 
jects (brown dwarfs and giant planets) can be pr edicted follow- 
ing th e non-gray theory of substellar objects by iBurrows et al.l 
( Il997h (B97 in the following), who have computed IR magni- 
tudes of substellar objects depending on age and mass. The up- 
per mass limit of planets is not yet defined. There are basically 
two suggestions: Either all substellar objects below the D burn- 
ing limit at around 13 Jupiter masses are refeiTed to as planets 
or all objects below the radial velocity brown dwarf desert (at 
20 to 50 Jupiter masses) are called planets, i.e., either all objects 
below 13 Jupiter masses or below 50 Jupiter masses. Calculation 
like those by B97 start with assumed initial conditions about the 
internal structure and temperature of the substellar objects and 
should, hence, be taken with care for very young objects up to 
a few Myrs (for a discussion see Baraffe et al. 2002). However, 
we use these calculations for our estimates due to the lack of 
better suited ones. According to B97's calculations, at an age 
of 100 Myrs and a distance of 300 pc, a companion at the 
Deuterium burning mass limit (~ 13 Mjup) would have an H- 
band magnitude of 22 mag - easily achieved by current opti- 
cal telescopes. We present in the following the first results of 
our search for substellar companions ar ound NSs. We note tha t 
part of our data have been analysed by iLo Curto et alJ (l2007h . 
However, they were interested in the spectra of the neutron stars 
as well as in hypothetical fallback disks around them. They nei- 
ther found a fallback disk nor was the IR data deep enough to 
constrain the spectra of the NSs. 



2. Sample selection 

While substellar companions may survive the supernova (SN), 
this case is probably rare, thus formation after the SN likely 
for most substellar companions. Hence, we use the NS age as 
good approximation of the age of possible substellar compan- 
ions. In case of SN survival the companion's age can be sig- 
nificantly underestimated. All NSs younger than 100 Myrs and 
closer than 300 pc were chosen as the initial sample in 2000. The 
radio pulsars amongs t these neutr o n star s were selected from 
the pulsar catalog of ' Taylor et~an (Il993[) . They were thought 
to be ordinary radio pulsars, thus still young as radio emission 
from non-recycled NSs is expected to last only at most 10 Myrs 
jLorimer & Kramej2004l) . To that sample, the radio-quiet X-ray 
thermal isolated neutron stars (XTINSs) were added. As noted 
above, the XTINSs are especially interesting to constrain the 
equation of state due to their blackbody-like spectra resulting 
in the best radius estimation currently available. They also rep- 
resent nearly half of the local y oung neutron star population (for 
reviews on these NSs, see, e.g., lKaplanll2008l: lHaberlll2004 . No 
close co-moving companions down to substellar masses have 
been identified in deep optical observations and proper mo- 
tion studies (e.g..lMignani et al.l2003l: I Walter & Matthewsll997t 
iMotch et al.ll200l 



The ages of the radio pulsars in our sample are estimated 
roughly from the characteristic spin-down age. However, the 
characteristic age is only a very coarse estimation of the true neu- 
tron star age. For the well known Crab pul sar the spin-down age 
exceeds the true age by about 25% (e.g., Manchester & Taylor! 
Il977l) . Initially assumed to be younger than 100 Myrs, some 
of our selected radio pulsars seem to be older from newer ra- 
dio timing data by Hobbs et al. (2004), so that possible sub- 
stellar companions would be fainter (see individual discus- 
sion in Section |5]l. The additional radio-quiet XTINSs found 
as soft ROSAT sources are most probably also young NSs 
(< 100 Myrs): Their proper motions indicate likely origins in 
nearby young star-forming clouds. The XTINSs have high mag- 
netic fields, an indicator of youth, and inferred spin-do wn ages 
from X-ray pulsations are less than 10 Myr (see, e.g., KaplanI 
2008). Thus, the XTINSs are probably several Myrs old, with 
cooling surfaces which emit the X-ray emission. For all seven 
radio-quiet X-ray thermal NSs, it appears unlikely that they are 
old with their X-ray emission due to accretion from the inter- 
stellar medium (ISM). First, they have constant X-ray emission 
over years wi th one notable exception where precession of the 
NS is likely ( Haberl et al.ll2006h . iNeuhauser & Trumped (Il999h 
quoted also statistical reasons for excluding accretion from the 
ISM. Furthermore, three radio-quiet X-ray thermal NSs have 
known proper motions, which are so high that they imply too 
low Bondi-Hoyle accretion efficiencies to describe these stars as 
old (10^) NSs re-heated by accretion from the ISM. 

For the radio pulsars, rough distance estimates are avail- 
able from the dispersion measure, which indicates the elec- 
tron column den sity towards the pulsar jTaylor et al.l 119931 : 
M anchester et al.l [2005). For the radio-quiet X-ray thermal 
NSs, a parallax measurement is available only in two cases, 
na mely RX J0720.4-3125 a nd RX J 1856.5-3754. For the lat- 
ter. I Walter & Lattimeii (|200^ ) determined a distance of ~ 1 17 pc 
with the Hubble Space Telescope. Recently Ka plan etaP (l2007h 
reported an improved distance estimation w ith RX J 1856.5-3754 
being at 167+|^ pc. For RX J0720.4-3125 lKaplan et all (l2007h 
reported a parallactic distance of 36 1;^^^^^ pc. The measured col- 
umn densities from the X-ray spectra of the other XTINSs are 
similar to those of RX J0720.4-3125 and RX J1856.5-3754, all 
smaller than the Galactic values (towards the respective direction 
of sight). One can assume the same order of magnitude for their 
distances as for RX J1856.5-3754and RX J0720.4-3125. Taking 
into account absorption by the inhomogenously distributed ISM, 
the estimated distances for tw o further radio-quie t X-ray thermal 
NSs support this assumption dPosselt et aLlbOOTl) . Here, we take 
300 pc as lower distance limit for the other 5 XTINSs. 
In the following we will present the observations and results for 
nine neutron stars, four of them observed at two epochs. All ob- 
jects, their known or assumed distances, and their proper mo- 
tions are listed in Tab. [T] 

Given the age and distance d of the NSs, the expected i/-band 
magnitude to find brown dwarfs and massive gas planets can be 
derived applying the results by B97. Bolometric corrections in 
the ^T-band h ave a reported range fro m a; 2 (T5 dwarf) to w 3.5 
(L5 dwarf) (Golimows ki et al. l l2004t) . We use in the following 
always a bolometric correction, BC = BCh - BCk = 3.5 mag. 



3. Observation summary and data reduction 

3.1. First epoch observations 

Nine neutron stars were observed in //-band utilising the Very 
Large Telescope (VLT) equipped with Infrared Spectrometer 
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And Array Camera (ISAAC) in service mode. The objects and 
the details of observations are listed in TablesfTland l A. 1 | in Sec.lAl 
of the Online Appendix. 

The VLT observation durations ranged from two to 
three hours. The seeing was usually better than 0.8" and 
the airmass in most cases below 1.5. The jitter-Mode was 
used with individual exposures of 12 sec and a jitter box 
width of 20". The data have been reduced with E SO's 
Eclipse (version 4.8.1 and 5.0.) dPevillardl [T999ll2001h tak- 
ing into account the corresponding darks, sky flats and bad 
pixel maps obtained during ESO's regular calibration plan. 

Table 1. The observed neutron stars 

Listed are for each object; the distance D, the proper motion 
in Right ascension {fj.„) and Declination (fig) if known (n.k. 
stands for not known). 



Object 


D 


Ha 






[pc] 


[mas yr~'] 


[mas yr"'] 


PSR J0108-1431 


130^ 


< 26^ 


< 78^ 


RX J0420.0-5022 




< 138*' 


< 138'' 


Geminga 


250+if c 


142.2+1.2'= 


107.4+1.2'= 


RX J0720.4-3125 


3611?'^ 


-93 + 1.2'i 


52.8 + 1.3'^ 


RX J0806.4-4123 




< 16^ 


< 16^ 


RX J1856.6-3754 


167+!«d 


326.7 ± 0.8 


-59.1+0.7*= 


PSR J1932+1059 


361X^"f 


94.03 ±0.14f 


43.37 + 0.29^ 


PSR J2124-3358^ 


270 ±20^ 


-14.4 + 0.8'^ 


-50 ± 2'' 


RBS 1774"^ 




n.k. 


n.k. 


^ IRXS J214303 


7+065419 






A PSR J2124-3358 turned out 


to be too old for 


our substellar 



companion search. See Online Appendix Sec.|C]for a NIR 
Umit of th is NS star 

^ i Mignani e t al. (200|), ^ 2cr Umit of iMotch et al.1 (l2006h . 

*= Fahertv et al. (2007), '''Kaplan et al.' ('2007'), 
Walter & Lattimer (2002),'^ Chatterjee et al] ( |2004 . 

g lOaensler et all (120021) . iHotan et al.1 (l2006l) 

^^'^ notes an assumed distance, which is a lower limit as 
indicated by estimatio ns from X-ray mea sured hydro- 
gen column density (e.g. lPosselt et al.l2007h . see also text. 

Individual images with poor seeing or high airmass have been 
excluded. For comparison also ESO's ISAAC pipeline has been 
applied. Multiple observations of the same object have been 
reduced individuall y. The results have then been combined using 
IRAF (Ijodv 19861 Il993h if observing conditions were similar. 
We fitted the absolute astrometry for the first ep och observation 
with the 2MASS point source catalogue (PSC, ISkrutskie et al.l 
l2006l) us ing the Graphical A stronomy and Image Analysis Tool 
(GAIA) (lDraperetal.ll2007h . In case of many relatively bright 
2MASS sources in flie fie ld (e.g. RX J1856.6-3 754), a positional 
accuracy of 300 mas (3o-. ISkrutskie et al.ll2006h can be achieved. 
In some fields only few 2MASS sources are appar ent. In these 
cases we used the USNO B 1 catalogue for which iMonet et al.l 
(1200 31) stated an astrometric uncertainty of 200 mas for absolute 
astrometry. However, this optical catalogue is based on much 
older observations than 2MASS and its entries must not nec- 
essarily coincide with apparently matching NIR detections in 
our observations. We use absolute astrometry only to derive 
NIR limits for the neutron star position taking an appropriate 
error circle to consider astrometric calibration effects as well as 
proper motions of the neutron stars. Neutron star positions are 
preferably either a VLBI position (e.g. PSR J1932+1059) or the 



position of the optical counterpart (e.g. RX J1856. 6-3754), if 
those are not known the (XMM-Newton) X-ray position (e.g. 
RBS 1774) is taken, references are corresponding to those of the 
proper motions in table[T]or those of the individual introductions 
in Sec. 15.1.1] to 13. II The NS's positional uncertainties have no 
impact on our search for substellar companions; in Sec. 15.31 
appropriate sky regions are considered for a possible NIR 
counterpart to a NS. 

No standard stars were explicitely observed for this project 
since only rough photometry was needed. We applied for photo- 
metric calibration the zeropoints obtained by ESO during their 
regular ISAAC calibration plan as starting point, and calibrated 
further with 2MASS photometry by using - preferably high 
quality - 2MASS point sources in the field (for individual de- 
tails see Sec. Is). 



3.2. Second epoch observations 

All the second epoch observations (see Tab. lA.ll ) were also car- 
ried out in the //-band in a similar jitter mode as the first epoch 
observations. The second epoch observations of RX J 1856. 6- 
3754 and PSR J1932+1059 in May 2006 were moved by ESO 
to the 3.5 m New Technology Telescope (NTT) equipped with 
the Son of ISAAC instrument (SOFI) (instead of VLT ISAAC). 
The two fields were observed with S OFI's Small field objective 
for around three to four hours during each of the four nights. 
The seeing was usually better than 1", and during the second 
and third night most of the time even better than 0.8", and in the 
last night very variable. The observations were reduced by ap- 
plying the S OF I pipeline procedures at the telescope. We cross- 
checked the result by reducing individual runs also with Eclipse. 
The further data reduction procedure was the same as for the 
first epoch observations applying IRAF and GAIA. The second 
epoch S OF I observations were up to 1 mag less deep than the 
ISAAC observations. 

In May 2007 VLT NAos COnica {NACO) observations were 
done for the inner field around PSR J1932+1059 using the 
coronographic mode for a bright star close to the pulsar NACO 
is an adaptive optical (AO) imager and spectrometer and has 
been operated with the S54 camera. The Visual dichronic ele- 
ment and wavefront sensor (0.45 -1.0 pm) and the opaque 1 .4" 
coronograph were used. A jitter box of 20", coronographic in- 
dividual exposures of 60 s, sky individual exposures of 0.3454 s, 
4 "object" offset positions and 6 "sky" offset positions were fur- 
ther parameters of this NACO run. Observing conditions for the 
roughly one hour observation were very good with seeing better 
than 0.6 at an airmass of around 1.3 and a reached strehl ratio of 
around 17% . The NACO observations have been reduced ap- 
plying the jitter engine of Eclipse (for more details, see SecfAl 
of the Online Appendix). 

Second epoch observations of RX J0720.4-3125 and Geminga 
were obtained with VLT ISAAC in December 2006 and 
January/February 2007 respectively (for more details see 
Tab. IA.lt . The data reduction procedure was the same as for 
the first epoch ISAAC observations applying Eclipse, ISAAC 
pipeline, IRAF and GAIA. 

The NACO field of view is 56" x 56" and therefore smaller than 
those of ISAAC and SOFI with 154" x 154". The (nominal) 
pixel sizes are 0.054" per pixel for NACO, 0.148" per pixel for 
ISAAC, and 0.144" per pixel for S OF I. 
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4. Method 

From the known proper motions of the neutron stars, it is possi- 
ble to calculate the expected shifts in pixel coordinates. First and 
second epoch observations have to be compared regarding object 
detections, magnitudes and positions. Due to the chosen time 
span it should be possible to distinguish clearly the expected 
shifts for co-moving objects from position errors. We tested sev- 
eral source detection algorithms, e.g. the Source Extractor by 
^ertin & Arnouts ( 1996), to determine the best relative astrom- 
etry. We refer to Sec. |B]for more details. Here we concentrate 
on our fi nal approach for the rel ative astrometry - the StarPinder 
-code by iDiolaiti et al.l (l2000bl) . an IDL-based point source de- 
tection algorithm for crowded fields which starts by determin- 
ing the point spread function (PSF) from the science image for 
further fitting. The pe r forman ce of StarFinder has been tested, 
e.g., by IDiolaiti et all (l2000al) and IChristou et all (12004) . The 
latter found, for example, that in fields with overlapping stellar 
PSFs StarFinder was capable to produce consistent astrometry 
of the order of one tenth of the PSF or b etter Photometric er- 
rors are usually bette r than around 0. 1 mag ( IDiolaiti et al.l2000at 
IChristou et al.ll20()4l) . For very faint sources StarFinder is sensi- 
tive to the structures in the determined model PSF. We carefully 
selected the reference sources for the model PSF and noise es- 
timation. Only objects detected with a SNR of at least three in 
both epochs - ISAAC, SOFI and NACO observations - were 
considered for the search for co-moving objects. The corre- 
lation factor for the StarFinder detections were usually better 
than 0.7. All sources have been manually inspected and appar- 
ent galaxies or artefacts have been excluded. If we compared 
sources detected by two different instruments, then the {S OF I, 
NACO) pixel coordinates were converted to coordinates within 
the ISAAC reference frame. Source lists of two epochs were 
then correlated and we derived the positional shifts as presented 
in Section|5] In case of two epochs of ISAAC observations, we 
minimized systematic instrument-related errors like field distor- 
tions by a similar observation setup. In case of different instru- 
ments/telescopes, the additional systematic errors showed up in 
a broader scattering of the computed positional shifts. In case 
of the NACO observations this extra error is negligible in the 
chosen ISAAC reference frame. 

If no second epoch observation was obtained, e.g., due to 
still unknown proper motion, we investigated the NIR sources 
in a 5000 AU circle around the NS position for their optical- 
NIR colors to identify potentially interesting objects for follow- 
up observations. We refer to Sec. IB. 21 for more details on the 
procedure. 

At the positions of the NSs, we checked also for faint NIR- 
objects which could be unexpected possible counterparts of the 
NSs themselves. We determined an upper //-band limit for the 
NS if there is no NIR object within the positional error by tak- 
ing the faintest, closeby object having a SNR larger than 3. We 
checked, whether we could reach deeper NIR-limits at the po- 
sition of the NSs by combining first and second epoch observa- 
tions. 



5. Results 

5.1. Pixel shifts between two epochs of fast moving neutron 
stars 

All faint NIR sources near the known NS position were analysed 
as outlined in Sec. |4| Co-moving companion candidates, found 
by comparison of first and second epoch observations, would 



have a high probability to be true companions, because NSs have 
very large proper motions, and it would be unlikely that an unre- 
lated, distant, faint background object would share such a large 
proper motion. 

Table |2| lists some parameters of the analysis. Additional 
notes on these and the results can be found in Sec. |5.l.l1 to |5.l.4| 



5.1.1. Geminga 

Geminga is the proto-type of the y-ray pulsars (see e.g. 
Big nami & Caraveo 1996). Th e distance of Geminga is still 
under discussion. FFahertv et al.l (|2007) recently reported a newly 
estimated parallactic distance of 250^g2" P'^ comparison with 



the older value of 157+;^^ pc by Ca raveo et all ( Il996lll998l). We 
adopt in the following the most recent distance by Fahertv et al.l 
(l2007h . Geminga has a spin-down age of around 340000 
years . The lowest age c o nsider ed by B97 (or similar papers, 
e.g., iBaraff^e et al.l I1998L |2002|) is 1 Myr We consider 383 
sources in the whole ISAAC field of Geminga, most of them 
fainter than 18.2 mag. The pixel shifts of the objects, common 
in both epochs, are plotted in Fig. (Tj There is no apparent 
co-moving source within w 10000 AU around Geminga down 
to // = 22.4 mag. For a 1 Myr age (see above) and 250 pc this 
corresponds to exclusion of all young companions having more 
than 12 Jupiter masses following the calculations by B97. 



5.1.2. RX J0720.4-3125 

RX J0720.4-3125 is the second brightest of the X-ray thermal 
isolated neutron stars. We assume in the following a similar 
age - 1 Myr - as for RX J 1856.6-3754. The distance to RX 
J0720.4-3125 has still a large uncertainty. iKaplan et all (l2007l) 
reported a parallactic distance from optical observations to 
be 361^Jg^ pc. The hydrogen column density fitted in the 
X-ray spectrum indicates a distance of 240 to 270 pc taking 
into acco unt the local inhom ogenously distributed interstellar 
medium (iPosselt et al.ll2007l) . In the second epoch, we had to 
reduce the correlation coeflicient constraint to 0.65 to re-identify 
all first epoch detections. 14 sources (out of the considered 
581) in the second epoch have correlation coefficients between 
0.65 and 0.7. The pixel shifts of the objects common in both 
epochs are plotted in Fig. [1] There is no apparent co-moving 
source within « 20000 AU around RX J0720.4-3125 down 
to // = 22.9 mag. For a 1 Myr age and 361 pc distance, this 
corresponds to exclusion of all young companions having more 
than 15 Jupiter masses following the calculations by B97. 



5.1.3. RXJ1 856.6-3754 

RX J 1856.6-3754 is the brightest of the X-ray thermal neu- 
tron stars. Its HS T- d eterm ined parallax of 117 + 12 pc 



by IWalter & Lattimed (|2002|) . was revised to 142+;^^ pc by 



iKaplan e t al.' ("2002") and recently further revised to preliminary 
167;|;[g pc (Kaplan et al. 2007). We adopt in the following this 
most recent value. The age of RX J 1856.6-3754 is assumed to be 
around 1 Myr gi ven the proper motion and the likely birth place 
in Upp Sco (e.g. I Walter & Lattimeij|20()2l) . There are 525 com- 
mon objects having a SNR larger than three in both, the ISAAC 
and S OFI observations, 448 of them fainter than H - 17.3 mag, 
and two with H > 21.4 mag. In Fig. [T]we plot the pixel shifts 
between the two epochs as well as the pixel shift of the neutron 
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Table 2. For objects with second epoch observations the following parameters are listed in this table: the brightest //-band mag- 
nitudes expected for a brown dwarf companion (mgo) or a gas planet of 13 Jupiter mass (mi3y„p) taking into account the NS's 
respective distance and age (see text and e.g. table [T]). Furthermore, the reached magnitude limit( m/„„) in both epochs, the time 
difference between first and second epoch, the covered projected physical separations iXproj), and the pixel scale of the reference 
ISAAC image (refX, refY) are listed. Reference image axes are roughly aligned with RA (X) and DEC (Y) within 1 degree, however 
sometimes with opposite direction (indicated by '-' instead of '+'). 



Neutron Star 






m;„„ 


epoch diff. 




refX 


refY 


aligns 


aligny 




[mag] 


[mag] 


[mag] 


[yr] 


[AU] 


[mas] 


[mas] 


+/- 


+/- 


Geminga 


18.2 


22.3 


22.4 


3.05 


10000 


148.4 


148.4 


+ 
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Fig. 1. Pixel shifts between the first and second epoch NIR observations of Geminga, RX J0720.4-3125, RX J1856.6-3754 

(from left to right). 

Errors are 3cr-StarFinder positional errors. The grey positions marked with Geminga, RXJ0720, RXJ1856 show the positions a possible 
co-moving object would have. For more details regarding the ISAAC reference frame pixel scale, axes orientation and covered magnitudes see 
Tab.|2]as well as Sec. |5TT] to [5T3l 



star which would be the same for a co-moving object. We did not 
account for pixel scale errors or instrumental distortion errors. 
Both are present in the scatter distribution of the object's pixel 
shifts. We do not expect an additional positional eiTor larger than 
one fifth pixel (mean 3cr positional uncertainty) due to coordi- 
nate transformation (SOFI to ISAAC reference frame). Even 
accounting for this does not change the picture of Fig. [1]- there 
is no co-moving //-band source with SNR > 3 present in our ob- 
servations. Applying the calculations of B97 for an assumed age 
of 1 Myr and a distance of 167 pc, substellar companions heav- 
ier than 0.01 Mq or 10.5 Jupiter masses can be excluded around 
RXJ1856.6-3754. 

A short note about the two objects with exceptional pixel 
shifts in the right-hand panel of Fig. [T] The object with offset 
of about 1.8, 0.5 belongs to a triple ISAAC source. It has an 
ISAAC //-band magnitude of 20.6 mag, the other two sources 
have H = 19.6 mag and 21.3 mag and are located roughly 1.2" 
from each other. In the S OF I observations, these three sources 
are more bluiTed and StarFinder can detect only two of them 
with SNR> 3 and correlation coefficient > 0.7. The derived sec- 
ond epoch position of the object in question is slightly shifted to- 
wards the position of the not detected fainter object. There were 
around ten such cases of non-detected multiple sources in this 
particular field. Usually, the positions derived by StarFinder are 
only marginally influenced. The pixel shift of around 2.5, 4.4 be- 



longs to an 21.0 mag object which is separated around 1.6" from 
a 16.5 mag object. By eye inspection this source is barely visible 
in the S OF I observation. StarFinder manages to still detect the 
faint source, however, the position seems not to be accurate for 
this high contrast. 

5.1.4. PSR J1 932+1 059 

PSR J 1932+ 1059 is a radio pulsar. New VLBI measu rements put 
this n eutron star at a distance of 361 ± 10 pc ( Chatterjee et al.l 
12004 . further away than previously thought. Therefore, only 
brown dwarfs could be detected as possible substellar compan- 
ions given its spin down age of 3.1-10'' years (Hobbs et al. 20041), 
and foflowing the models by B97. PSR J1932+1059 hes unfortu- 
nately in a crowded field with a very bright, // = 8.8 mag, object 
only X 2.5" away from the pulsar. The combined S OF I obser- 
vations of May 2006 were more blurred than the ISAAC ones 
due to the observing conditions (small visibility time window, 
not so good weather on the last night). The bright NIR object 
causes large "doughnut"-shaped artefacts (only in the S OF I im- 
age) due to the jittering process. While the StarFinder -Code is 
able to still detect reasonable well the NIR sources, we avoided 
most of the artefacts and the bright source by using the NACO 
coronographic observations from May 2007 instead. Note that 
the NACO S54 camera has a pixel scale of around 54 mas com- 
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Fig. 2. The A^ACO-field of PSR J1932+1059 with circles indicating common sources in NACO and ISAAC observations. The big 
boxes are pointing toward muhiple NACO detections which have been identified as only one source in the ISAAC image. Triangles 
with ' S ' mark positions of sources for which we used the S OF I observation for the common proper motion search. The right figure 
is a zoomed-in (roughly 14" x 10") variant. Sources not considered for the common proper motion search are marked with squares. 
They are detected only in the NACO image due to the coronographic mode, but not in the ISAAC observation. The (empty) position 
of the pulsar is indicated by the diamond shape with 1" radius. 



paring to around 148 mas for ISAAC. NACO StarPinder posi- 
tional errors are much smaller than positional errors for ISAAC. 
NIR sources in an area of around 10000 AU projected physi- 
cal separation of the pulsar (at an assumed distance of 361 pc) 
have been considered for further analysis. The area covered by 
NACO extends in the east and north up to about 3000 AU, and 
in the west and south up to roughly 6500 AU projected physical 
separations from the NS. In this area, we used the SOFI im- 
age only for four sources, which we could not detect well in the 
NACO image (wire construction shadow). Due to the corono- 
graphic NACO AO mode and smaller pixel scale there were 
more sources detected as in the ISAAC image. By checking the 
NIR sources found in both epochs, several single ISAAC sources 
turned out to be actually multiple NACO sources. StarPinder 
positional errors for the ISAAC sources are not conspicuously 
high in these cases indicating the limits of the StarPinder code. 
In case of a relatively bright main component in the multiple 
source, it is reasonable to assume that this specific source was 
actually detected in the ISAAC image. If the multiple NACO 
sources had nearly the same magnitudes we compared the flux- 
weighted average position of the individual NACO sources to 
the one ISAAC position in such cases. Furthermore, we looked 
at the pixel shift of each individual NACO source comparing to 
the one ISAAC source aiming to get a feeling how large pixel 
shifts due to unidentified multiple sources could be. All con- 
sidered common NACO-ISAAC sources are marked in Fig. |2] 
Close (<^ 4") to the bright star there are eleven NACO sources 
which were not detected with ISAAC at all. These sources have 
//-band magnitudes of around 18 to 21 mag. They are marked 
in Fig. |2]and have to be investigated for co-motion in another 
epoch. 

The pixel shifts of the 284 sources considered for the com- 
mon proper motion search in the SOFI and the ISAAC ob- 
servations are plotted in Fig. [3] Due to the observing condi- 
tions the errors in the SOFI observations are relatively large. 
There is a faint 21.2 mag object labeled 'XT in Fig. |3](left) 
whose 3cr error bars overlap with the position of an object which 



would share common proper motion with PSR J1932+1059. The 
spread in the pixel shifts in Fig. [3] especially the fragmenta- 
tion in two 'clouds' indicate systematic effects. The SOFI small 
field is known for strong distortions in some parts of the in- 
struments field of view. This may explain XI 's position in the 
pixel shift plot. To account for distortions we applied SCAMP 
(Bertin 2006) with a polynomial degree of five for both, ISAAC 
and S OF I images. As expected there are large S OF I distortions 
while for ISAAC they are small. We applied then S Warp Q by 
E. Bertin (Terapix project) for resampling and proceeded with 
StarPinder for a similar analysis as described above. The re- 
sulting pixel shift plot is shown in Fig. [3] (center). The sources 
are much more clustered in this figure. Sources with large pixel 
shifts are either indeed moving objects or they are close to bright 
stars with their detection parameters possibly slightly affected. 
The Xl-object and the pulsars proper motion do not coincide 
within their 3cr uncertainties, however only marginally. It is very 
likely that XI is a background object. If unexpectedly XI is a 
co-moving companion, this source would be at around 9200 AU 
projected physical separation from the pulsar, and applying the 
calculations by B97, it would have a mass of around 42 Jupiter 
masses for a distance of 361 pc and an age of 3.1 ■ 10^ years. 
There is no known object within 10" of XI in the SIMBAD- 
database. Unfortunately, XI was not in the field of view of our 
NACO observations. 

The pixel shifts of the 104 NACO sources with respect to 
the earlier ISAAC observations are shown in Fig. [3] The NACO 
sources 'I a' and 'I b' belong to one 17.2 mag ISAAC detec- 
tion as well as 'II a' - 'II c' to another 19.1 mag one. Here, 
it seems likely that the single ISAAC sources are actually a 
merger of the multiple NACO ones. The flux-weighted average 
NACO positions can be regarded as the putative ISAAC posi- 
tions if this instrument would have been used. The respective 
flux- weighted average pixel shifts 'I weigh' and 'II weigh' are 
situated within the bulk of the pixel shifts of the other sources. 



' http://terapix.iap.fr/rubrique.php?id_rubrique=49 
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Therefore, the ISAAC sources T and 'IF can be ruled out as 
co-moving objects. The large pixel shifts of I a-b, II a-c (each 
NACO source considered separately) illustrate that analysing de- 
tections by two instruments with different spatial resolution is of 
limited use given the high precision needed for the relative as- 
trometry. Thus, to sufficiently probe the co-motion of all individ- 
ual NACO sources another NACO epoch would be necessary. 
Overall, Fig. [3] shows for the common NACO-IS AAC sources 
that there is no pixel shift which would be expected from a co- 
moving object. However, another NACO observation would be 
needed to confirm this. 

The current state of our search around PSR J1932-H1054 is 
that we found no co-moving object down to // = 21 .5 mag com- 
paring the ISAAC with the SOFI observation, as well as the 
ISAAC with the NACO observation. At the given distance and 
age, there are probably no substellar companions down to 42 
Jupiter masses. However, this result has to be seen in the light of 
the described uncertainties due to the different instruments and 
the close bright star. Further investigations are needed for con- 
firmation. 
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Fig.S.FaintNIR-objects ai-oundRXS J214303.7+065419 which 
may be interesting for second epoch observations. The big circle 
would have a radius of 5000 AU at 300 pc. The marked faint 
objects have a detection SNR of at least three. 



5.2. Neutron stars with only first epocli observations 

In the following, the first epoch //-band observations in connec- 
tion with available optical data are discussed to judge about the 
necessity of second epochs to the ISAAC observations. 

PSR J01 08-1 431 is one of the closest radio pulsars accord- 
ing to its very low d ispersion measure (2.38 + 0.01 cm-^^ pc, 
iD'Amico et al.lll998l) . The faint p ulsar has a distance o f 130 pc 
and seems to be older (170 Mvrs. iMipnani et aLll2003h than we 
assumed initially. Possible brown dwarfs around PSR J0108- 
143 1 are expected to be fainter than H - 20.8 mag. The obtained 
photometric magnitudes for the objects of interest have large er- 
ror bars due to the availability of only two bright 2MASS sources 
for the zeropoint determination. Astrometric positions were re- 
fined using USNO Bl point source positions. Within 5000 AU 
there is only one source fainter than H - 20.8 mag having 
H = 21.0 mag. We obtained the archival VLT FORS 1 images in 
the V-band, reduced the data applying bias subtraction and flat- 
fielding and calibr ated the source V-ban d magnitudes with the 
tabulated values of iMignani et all (l2003h . The faint NIR-object 
is detected in the V - band with 25.6 mag. Therefore, the V -H- 
color is 4.6 mag. This color is typical for spectra l type M2 V 
or M3 V according to Kenvon & Hartmamil (Il995h . Even when 
taking an error of 1 mag into account the source is no candidate 
for a substellar companion around PSR J0108-1431. Thus, there 
are no substellar companion candidates for second epoch obser- 
vations from our VLT ISAAC image within 5000 AU around 
PSR J0108-1431. 

RX J042 0-5022 is one of the X-ray thermal isolated neutron 
stars (e.g. iHaberl et al.ll2004l) . Neither distance nor proper mo- 
tion were measured due to a still insecure optical counterpart. 
The hydrogen column density fitted in the X-ray spectrum indi- 
cates a distance of very roughly 300 pc. Assuming a similar age 
as that of the better known XTINS RX J 1856.6-3754, one would 
expect substellar companions to be fainter than H = 19 mag 
at this distance. Only two objects are present within 5000 AU 
around the neutron star position for an assumed distance of 
300 pc. There are optical archival observations with ESO's VLT 
equipped with FORS 1 (66.D-0128). The three B-band observa- 



tions were added and the photometry was calibrated with the 
USNO Bl catalogue due to missing archival calibrations. At 
the position of the two interesting sources there are no optical 
sources down to the w 26 mag limit of the FORS 1 B-band ob- 
servations. Thus, the NIR-sources could not be excluded as pos- 
sible substellar companions by their B-H color. Only the second 
epoch observation, as soon as the proper motion is constrained, 
can clarify possible companionship. 

RX J0806.4-4123 is one of the faintest of the seven XTINSs 
(e.g. Haberl et al. 2004). Neither proper motion nor distance 
were measured due to a still missing optical counterpart. The 
hydrogen column density fitted in the X-ray spectrum indicates 
a distance of 240 + 25 pc taking into accoun t the local imhomge- 
nously distributed ISM (Posselt et all 120071) . Possible substel- 
lar companions around RX J0806. 4-4123 should have //-band 
magnitudes fainter than 18.6 mag We find around 50 interesting 
NIR-sources within 5000 AU down to H - 22.0 mag having a 
minimal SNR of three. We observed the very crowded field in 
the V-band using ESO's NTT equipped with the Superb S eeing 
Imager -2{S US 12) in December 2004. The data were reduced 
applying dark field subtraction and flat fielding. Bias subtraction 
was neglected. Astrometry of the optical observation was ob- 
tained by fitting with the USNO Bl point source positions. The 
photometry was obtained by a standard field observation reduced 
in the same way as the science field. The derived limit for 3cr de- 
tections is V = 24.8 mag. A few NIR-sources could be excluded 
as possible substellar companions due to their low V - //-color. 
However, most of the interesting objects are not detected in the 
V-band. All the final 45 objects will have to be investigated in 
the second-epoch observation. 

RXS J21 4303.7+06541 9, also called RBS 1774, is a XTINS. 
Currently, we do not know of an optical counterpart down to 
r' - 25.7 mag jRea et aLll2007h . The most recent X-ray position 
by Chandra has an uncertainty radius of 0.6" ( Rea et alj|2007i) . 
Here, we assume as a lower distance limit 300 pc and an age sim- 
ilar to that of RX J 1856.6-3754. Substellar companions are then 
expected to be fainter than H - 19 mag. There are nine inter- 
esting objects within 5000 AU around the neutron star position. 
Due to the availability of only three relatively bright and clus- 
tered 2MASS point sources the photometric error of the sources 
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Fig. 3. Pixel shifts between the first and second epoch NIR observations of PSR J1932+1059 

Errors are 3cr- StarFinder positional errors. 'PSR1932' marks the position of an object which would co-move with PSR J 1932+ 1059, 
the shown error bars correspond to the 3cr uncertainty of the pulsars proper motion. The left plot is for a distortion-Mncorrected 
analysis of the ISAAC and SOFI images. The 21.2 mag object labeled 'XT is 25.6" north-east of the pulsar The middle plot 
shows the pixel shifts of the same objects, but for distortion-corrected ISAAC and SOFI images. The re-sampled ISAAC pixel 
scale is 148.4 mas for x, and 148.4 mas in y. The object XI and the pulsars proper motion do marginally not coincide within 
their 3cr uncertainties. Pixel shifts for common sources in the ISAAC and the NACO observation are shown in the right plot. As 
described in the text single ISAAC sources turned out to be multiple NACO sources. Pixel shifts of the two components of the 
ISAAC source T and the three components of the ISAAC source 'IF are labeled in the plot as well as the flux-weighted average 
positions of Ia,b and Ila-c. See Tab. |2]and text for more details. 




Fig. 4. Faint NIR-objects around RX J0420.0-5022 and RX J0806.4-4123 which may be interesting for second epoch observa- 
tions. The big circles have a radius of 5000 AU at the respective distances. The marked faint objects have a detection SNR of at least 
three. 



is around 1 mag, especially for faint objects. The 2MASS ob- 
tained astrometry was improved by using additional USNO Bl 
point source p ositions. Th e archi val ESO EMMI /?-band obser- 
vations by Zampieri et al.l ( 1200 ll) have a limit of R -23 mag. 
EMMI detected objects, some also present as USNO B 1 point 
sources, provided the R - H color for the NIR-sources. For five 
NIR-sources only the limit of R - 23 mag could be applied and 
thus no statement of a color in favour for substellar companion- 
ship is possible. The four sources with a known /?-band magni- 
tude can be excluded as possible substellar companions due to 
their R - H color, even considering a color magnitude error of 
2 mag. Two of the sources are elongated and are probably galax- 
ies. Whether the remaining three sources are co-moving objects 
can only be constrained with the knowledge of the still unknown 
proper motion. 



5.3. H-band limits at fhe neutron star positions 

In the following we list the //-band limits obtained from our 
observations. Since the main focus of our observations has been 
relative astrometry no standard stars have been observed and sky 
conditions were generally not photometric. Furthermore, indi- 
vidual runs from different days have been combined to detect 
faint sources. Zeropoints (and their errors) obtained during the 
standard calibration plan by ESO were taken to derive a rough 
photometric calibration. The zeropoints are from the same date 
as the individual runs (see Table lA.ll ) and were averaged accord- 
ing to the combination of the runs. 

The StarFinder code was applied for the photometric mea- 
surements resulting in relatively small - nominal - magnitude 
errors. Due to the observing conditions (usually clear or thin cir- 
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Table 3. //-band and mass limits obtained from the observations 

Listed are for each object: ///,„, indicates the faintest //-magnitude with a detection SNR of three, which one could measure at the position 
of the neutron star. The Icr error is a resuh of the StarFinder photometric error and the zeropoint error and an estimated/assumed error 
from comparison with 2MASS. See also text. ZP states the magnitude of the zeropoint, which were obtained in the frame of ESO's 
calibration plan. ""^ indicates that two epochs have been used. S E - 2MASS - limits lists the Source Extractor limits if using relative 
2MASS photometry. The assigned eiTor indicates poor or high number of availa ble, suited 2MASS P SC sources, but includes also 
detection routine uncertainties. For comparison, the next column lists the values by lLo Curto et alJ (l2007h . which they derived from our 
observations. They stated a general magnitude error of 0.1 to 0.2 mag. The last column notes the reached mass limit (from Sec. 15. lb . ■ ■ ■ 
notes that the found objects are interesting for a second epoch observation, - states the opposite. 



Object 


//,,„, ± icr 


ZP 


SE-2MASS - limits 


Lo Curto et al. f2007) 


Mass-limit 




[mag] 


[mag] 


[mag] 


[mag] 


[Jupiter masses] 


PSRJ0108-1431 


21.42 ±0.16 


24.66 


21.5 ± 1.5 






RX J0420.0-5022 


21.84 + 0.16 


24.61 


21.5 ±0.6 


21.7 




Geminga^ 


22.66 ±0.14 


24.66''2 


22.0 ± 0.4 




12 


RX J0720.4-3125 


23.07 ±0.11 


24.63'^ 


22.3 ± 0.4 


22.7 


15 


RX J0806.4-4123 


22.64 ±0.16 


24.59 


22.6 ± 0.4 


22.9 




RX J1856.6-3754 


21.54 ±0.24 


24.65''2 


22.1 ±0.4 


21.5 


11 


PSRJ1932+1059 


20.31 ±0.18 


23.49 






42 


RBS 1774 


22.04 ±0.16 


24.67 


21.6 ± 1.0 


22.1 





" Geminga has been detected in a deeper HST F160W observation with roughly 24.1 mag dShibanov et aU2()0^ 



rus) and the combinations of individual runs the error is larger. In 
Table [3] we report the //-band limits at the positions of the neu- 
tron stars as well as the determined minimal photometric errors. 
In case of a sufficient number (at least 5) of 2MASS sources with 
very good quality flags in all bands we compared our StarFinder 
magnitudes with the 2MASS values for the respective sources. 
A coiTection to our magnitudes was then applied as well as its 
additional magnitude error (standard deviation of the found dif- 
ferences). For fields without sufficient 2MASS source number, 
an assumed magnitude error of 0.1 mag (Icr) was considered 
additionally to the StarFinder and zeropoint error. This reflects 
the estimated uncertainties due to the sky conditions and the 
airmass correction for the respective fields. The average atmo- 
spheric extinction in the //-band at Paranal is 0.06 mag B We 
list also the Source Extractor magnitudes derived by only doing 
relative photometry with the available, preferably faint, 2MASS 
sources. We regard these magnitudes as much coarser since often 
only very few 2MASS sources are present (e.g. only two sources 
for the field of PSR J0108-1431) and not all have good quahty 
flags. Therefore we assigned a much higher magnitude eiTor to- 
wards these " SB 2MASS-limits". For some of our observations 
ILo Curto et a l. (2007) have already estimated //-band limits (by 
relative photometry with all 2MASS sources in the field) which 
are listed for comparison in Table [3] Except for RX J0720.4- 
3125, where we have a slightly deeper limit due to the second 
epoch observation, our values correspon d within the error bars t o 
the results for 4 XTINSs investigated byE o Curto et a D(|2007|). 

If two epoch observations are available, it would be interest- 
ing if the combined image gives a deeper limit. In case of PSR 
J1932+1059 the NACO observation improves the constraints in 
the //-band limit since the bright source close to the pulsar was 
behind the coronograph. In the NACO field sources are detected 
down to 21 .3 mag. Since photons from the bright star are present 
close to the position of the NS, this might hinder detection of 
very faint objects. Therefore we take as limit the H magnitude 
of the faintest source which has a similar angular separation from 
the bright source as the pulsar. Thus, the //-band limit for PSR 
J1932-I-1059 is 20.31 + 0.18 mag. 



In case of two ISAAC epochs (similar observing condi- 
tions) we combined both to test for a deeper limit for Geminga 
and RX J0720.4-3125. The Zeropoints are averaged for the two 
epochs and StarFinder measurements were refined with high- 
quality-2MASS point sources as described above. In case of RX 
J1856.6-3754 the SOFI image was transformed to the ISAAC 
reference frame and both images averaged. The StarFinder mag- 
nitudes were first calibrated roughly with the ISAAC zeropoint 
and then refined with 2MASS as well. In Table [3] we list the 
results of the combined observations. In Figure |6] the images 
around the NS's position are shown, in case of second epochs 
the combined image is chosen. The field of PSR J1932+1059 
is shown in Fig.|2] For RX J0806. 4-4123 one might suspect a 
very faint source within the r = 1.5" circle. However, this is not 
an 3cr detection (3cr //-band limit is 22.64 + 0.16 mag). Deeper 
observations would be of interest here although the field of RX 
J0806. 4-4123 is very crowded, and the chance of background 
objects is high. 



6. Discussion 

We analysed multiple-epoch //-band observations of four 
isolated neutron stars and found no co-moving object in our 
data. The derived limits in the frame of the B97 calculations are 
12, 15, 11, 42 Jupiter masses for Geminga, RX J0720.4-3125, 
RX J1856.6-3754, and PSR J1932-H1059 assuming ages of 1, 1, 
1, 3.1 Myrs and distances of 250, 361, 167, 361 pc, respectively. 
The small sample of investigated objects allows only crude 
statistical c onclusions. A new query at the most recent pulsar 
data base (" Manchester et al.l 120051 Fh for pulsars with ages 
below 100 Myrs and distances below 400 pc, results in ten NSs 
(including Geminga, and the XTINSs RX J0720.4-3125 and 
RX J 1856.6-3754). Listed among them is the double-pulsar 
PSR J0737-3039. However, we exclude this syst em since first, 
the distance stated in the corresponding reference ("Kram er et al.l 
2006) ranges from 200 to 1000 pc, and second, this would imply 
a merging of statistics for double NS systems with isolated NSs. 
We consider two XTINSs additional to RX J0720.4-3125 and 



' www.eso.org/instruments/isaac/tools/imaging_standards.html#BB 



" http://www.atnf.csiro.au/research/pulsar/psrcat 
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Fig. 6. Neutron star detection limits in the //-band 

Moving left to right, top to bottom the //-band fields of the following neutron stars are shown: first row: PSR J0108-1431, RX J0420.0-5022, 
Geminga, RX J0720.4-3125; second row: RX J0806.4-4123, RX J1856.6-3754, RBS 1774. Except one each circle has a radius of 1.5". The 
radius for the circle around the position of RX J1856.6-3754 has a radius of 2". The box sizes are roughly 16" x 13". See also text. 



RX J 1856.6-3754 as the upper limit for the number of XTINSs 
within 400 pc. This upper limit is chosen because first, already 
the second brightest XTINS, RX J0720.4-3125, has a distance 
of 361 pc, and second, we take into account the indications from 
the absorbing hydrogen column density and ISM models fo r 
the distances of the fainter XTINSs (e.g. iPosselt et al.ll2007l) . 
Thus, we consider 1 1 NSs within 400 pc as young (< 100 Myrs) 
enough to show co-moving companions with our technique. We 
note that we do not include the recently found Rotational Radio 
Transients, since the currently-ide ntified sources are thou ght to 
be at distances between 2 to 5 kpc (iMcLaughlin et alJ2OO60 . The 
probability, P, to get no detections in a subsample out of 1 1 NSs 
having an (unkown) likelihood, Lcom, for substellar companions, 
is equivalent to an hypergeometric distribution probability. P 
can be calculated according to the conventional equations of 
the hypergeometric distribution for different likelihoods Lcom- 
Given 1 1 sources, a zero detection rate of brown dwarfs among 
four NSs results in P =2%, 21 %, or 64%, if Lcom is 55%, 
27%, or 9% , respectively. For planets heavier than 15 Jupiter 
masses, the zero detection rate of substellar companions among 
three NSs results in P = 6 %, 34 %, or 73 %, if La,,,, is 55 %, 
27 %, or 9 % , respectively. It is likely that the young population 
of isolated NSs has a similar substellar companion fraction 
as older isolated NSs. Strictly spoken, the numbers presented 
above cover only the case of formation of substellar companions 
after the supernova since we approximated the oldest age of 
possible substellar companions with the NS age. However, the 
SN survivor's rate is likely to be even smaller. We conclude 
that brown dwarfs and giant gas planets heavier than 15 Jupiter 
masses are rare among isolated NSs. Possible reasons are: 1. 
These substellar companions did not survive the destructive 
past of the NS. 2. They are too massive to be formed from a 
supernova fallback disk. 3. Massive giant planets (> 10 MjupHer) 
are rare, at least among the well-studied low-mass stars. Less 
than 5% of the known planets are such objects (March 2008, 
SchneideiQ). 4. Only few massive stellar progenitors host giant 



' The Extrasolar Planets Encyclopaedia: |http://exoplanet.ei]71 



planets. There are indications that the rate of substellar com- 
panions decreases with ste llar mass. For examp le, considering 
low-mass (< 5 Mq) stars, [Kennedy & Kenvonr(l2008i) found a 
peak at stellar host masses ~ 3 M© for the probabiUty of having 
at least one gas giant around them. While from massive star 
formation sites one knows disks around protostars with masses 
less than ~ 20 M© (Cesaroni et al. 2007), these disk are expected 
to have a shorter lifetime than those around low-mass stars 
and this may reduce the formation rate of massive substellar 
companions. There is currently no planet known around a star 
more massive then 4.5 M© (March 2008, Schneider), which 
may be, however, partly due to instrumental limitations (e.g. 
brightness contrast). 5. Substellar companions may be fainter 
than calculated in B97. Comparing to other theoretical models 
based on non -gray atmospheres, the B97 effective temperatures 
are hotter ( Barafife et al.l l20(j2l) . Thus, actual //-magnitudes 
might be lower than calculated by B97. 

Seeing our result in the light of the searches for substellar 
companions around white dwarfs it is interesting that substellar 
companions seem to be sparse among white dwarfs too: until 
now only three L-type compa nions amo ng well over 1000 
white dwarf targets are known (Farihi et al.li2008. and references 
therein) and one s uspected planet (~ 2 M jupHer) around a pul- 
sating white dwarf (iMullallv et al .120081) has yet to be confirmed. 

A short note on RX J0720.4-3125: iHaberl et al.l (|2006|) 
suspected free precession to be responsible for the ~ 7 years 
variations in the pulsed X-r ay spectra. Note, that thi s is still 
a subject under discussion (Ivan Kerkwiik et al] l2007h . In the 
frame of our work it is interesting to note, why forced pre- 
cession, e.g., by a s ubstel lar companion, is not possible. For 
example, ILiu etld] (120071) discussed forced precession by a 
planet for PSR B1828-11. In t he case of RX J0720.4-3125, we 
can apply the same formula as iLiu et al. I (|2007) to approximate 
the orbital radius which would be needed for planets lighter 
than 15 Jupiter masses so that they could cause the precession 
observed. The largest orbital radius, such planets could have, 
would be 1 X lO"*" AU (for zero inclination, NS distortion of 
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4 X 10 ). However, this would be below the smallest Roche 
limit for a (normal) NS (radius of 10km, density of 10'^ kg m"-^) 
- both for planets with densities similar to Earth, resulting in 
a Roche limit of 2 x lO^^* AU; and for planets with densities 
similar to Jupiter, resulting in a Roche limit of 3 x lO^^* AU. 
Thus, even if there may be planets lighter than 15 Jupiter 
masses around RX J0720.4-3125, they cannot cause a ~ 7 years 
precession of the NS. 

In agreement with expectations no neutron star was detected 
in the //-band. The NIR limits are only weak constraints on the 
overall NS spectra given the known deep optical observations. 
This result is in agreeme nt with that of the inde pendent investi- 
gation for the XTINSs bv lLo Curto et all (|2007|) . 
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Appendix A: Observations 

The VLT NAos COnica {NACO) coronographic observations 
of PSR J1932+1059 were done with the S54 camera and the 
opaque 1.4" coronograph. As a reference star for the AO cor- 
rection the bright (V = 13.4) star GSC2 N0232103780 was se- 
lected, located 5" away from PSR J1932+1059. A jitter box of 
20", coronographic individual exposures of 60 s, sky individual 
exposures of 0.3454 s, 4 "object" offset positions and 6 "sky" 
offset positions were further parameters of this NACO run. The 
jitter engine of Eclipse was applied for data reduction consid- 
ering the corresponding dark and lamp flat, 28 individually se- 
lected objects for the image stacking algorithm and separate ob- 
ject and sky files. 

Due to the relatively small number of images, the leaking-in 
photons of the bright object outside the actual field of view, and 
the coronograph wire construction, the resulting image (Fig. |2]l 
shows several faint holes and artefacts which may result in false 
source detections. We inspected each found source and excluded 
such artefacts from further analysis. Since we are not interested 
in highly accurate photometry but in (relative) astrometry the 
image quality serves its intended purpose. 



Appendix B: Further details on methods 

B.1. Other source detection algorithms 

For the relative astrometry, we considered two d ifferent ap- 
proaches. First, we applied GAIA dPraper et alj|2007 ) for object 
detection provided that a signal-to-noise ratio (SNR) of at least 
three is reached. GAIA makes use of the Source Extractor by 
iBertin & ArnoutsI (11996) . In this program, which was designed 
with barycenter-detection of galaxies in mind, one can choose 
several detection parameters like the used filters (e.g. Gaussian), 
or how to handle overlapping sources. We tested several sets of 
parameters and applied individually the apparently best suited 
one for each observational field. We then estimated the //-band 
magnitude by calibrating with the Source Extractor (preferably 
faint) 2MASS PSC sources. The achieved SNR for a 22 mag ob- 
ject was around three, magnitude eiTors were around 0.4 mag. 
The positional accuracy for faint sources stated by the Source 
Extractor was often surprisingly good (better than one tenth of 
pixel for faint sources) which seems not realistic. We checked 
the positional accuracy reached by 2D-Gaussfits using IDL. For 
the interesting - usually faint - objects, three times the mean po- 
sitional uncertainty was 3cr ^ 80 mas, but could also reach ex- 
ceptionally 230 mas for the faintest. The pixel scale of ISAAC 
was » 147 mas. Brighter objects had usually better determined 
positions. While the 2D-Gaussfits were useful to illustrate the 
possible maximum positional eiTors, they failed for overlapping 
sources. Therefore, we applied as fi nal approach f o r the re lative 
astrometry the StarFinder -code by ' Diolaiti et al.l (l2000b ) with 
its known accuracy as noted in Sec.|4] For bright //-band objects 
the 3cr eiTor of the pixel shift is around one hundredth of a pixel, 
while for the faintest object it reaches one pixel. We compared 
the Source Extractor / GAIA source list with the StarFinder 
source list for a performance check. StarFinder was usually able 
to find at least as many or more sources as the Source Extractor 
of GAIA, and it only failed for galaxies which are excluded any- 
way from further analysis. 



B.2. Procedure in case of first epoch only observations 

If no second epoch observation was obtained, e.g., due to still un- 
known proper motion, we used the following procedure to iden- 
tify at least potentially interesting objects for follow-up obser- 
vations. Orbital separations of at least 5000 AU are inspected 
for NIR sources fainter than the magnitude obtained by the 
luminosity-magnitude relation and considering B97 at the brown 
dwarf upper mass limit. NS distances listed in Tab. [T] were used 
to determine the 5000 AU circle around the NS position. The 
relatively wide separations take into account recently found dis- 
tant substellar companions at roug hly 4000 AU around m ain 
sequence stars (see e.g. table 3 in iBurgasser et al.l (l2005h for 
brown dwarf companions) and possible widening of orbits due 
to the mass loss of the primary during the supernova according 
to Kepler's laws. 

The faint sources detected by GAIA were cross-correlated 
to sources in available optical observations obtaining optical- 
NlR-colors. This was usually done with the help of archival or 
kindly provided optical observations. Many first-epoch-objects 
could be rejected as possible substellar companion candidates 
due to these colors being typical for normal stars according to 
Kenvon & Hartmann ( 1995 ). The give n optical and optical-NlR- 
colors bv iKenvon & HartmannI {1995) go down to spectral type 
M6. As an example, one would expect the R - H color to be 
lower than 5 mag for normal stars down to M6 V, while L and T 
dwarfs usually have higher values o f around 6 mag. This can be 
seen in Fig. 1 of Kirkpatric k et alj (fT99 9) shown for the R - K 
color (the H - K colors are around 0.6 mag fo r most L and T 
dwarfs, see e.g. Fig. 14 of lKirkpatrick et aOl99^) . However, we 
could not exclude ultra-cool M-dwarfs (M8) at this stage. 

We note the use of a conventional magnitude error of usually 
1 mag up to exceptionally 2 mag for the optical-NlR-colors 
considering variations in the magnitudes of the reference optical 
observations as well as errors in the magnitude estimation of the 
NlR-objects. By this, all potentially interesting sources should 
be included. 
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Table A.l. List of our //-band observations 

The night parameters of each observation are averaged. Note that the DIMM seeing is the standard ESO V-band value. The Near 
Infrared seeing is usually better by at least 0.1 arcsec. In the column "Included ?" it is stated whether we took this observation run 
into account or did not use it, e.g., because of known technical problems or bad weather during the observations. Furthermore, it is 
noted by "Y-mod" if individual files were excluded from the observational run, e.g. due to exceptional bad seeing. For completeness 
we note that five neutron star fields were akeady observed in December 2000 using ESO's NTT equipped with the S OF I instrument 
in the //-band (66.D-0135). However, comparing to the later VLT ISAAC observations, the SOFI images were generally less deep 
and suffered from worse seeing conditions. We do not list the NTT S OF I observations from 2000 here. 



Object 


Instrument 


Exposure 


Nights 


Airmass 


DIMM Seemg 


Sky 


Included ? 






[s] 


yyyy-mm-dd 




[arcsec] 




Yes / No 


PSR J0108-1431 


VLT ISAAC 


1200 


2003-06-23 


1.1 


0.7 


photometric 


Y-mod 




VLT ISAAC 


1380 


2003-08-07 


1.0 


0.8 


photometric 


Y-mod 


RX J0420. 0-5022 


VLT ISAAC 


1980 


2004-01-10 


1.3 


0.6 


clouds 


N 




VLT ISAAC 


1980 


2004-01-13 


1.7 


0.9 


photometric 


Y-mod 




VLT IS AAC 


1980 


2003-1 1-13 


1.4 


0.8 


clear 


Y 


Geminga 


"\ 7T T C A A y — ^ 

VLT ISAAC 


1980 


2004-01-02 


1.6 


0.9 


photometric 


Y 




"\ 7T T C A A ' — ^ 

VLT ISAAC 


900 


2003-12-30 


1.4 


0.6 


cirrus 


Y 






LyoU 


zuuo- 1 i-i J 


1 A 


U.J 


clear 


V 

1 


RX J0720.4-3125 


VLT ISAAC 


1980 


2004-01-08 


1.0 


0.5 


cirrus 


Y 




VLT ISAAC 


1980 


2004-01-11 


1.0 


0.7 


clear 


Y 




VLT ISAAC 


1980 


2003-12-04 


1.0 


0.8 


photometric 


Y 


RX J0806.4-4123 


VLT ISAAC 


1980 


2004-12-21 


1.1 


0.7 


photometric 


Y 




VLT ISAAC 


1980 


2004-12-21 


1.1 


0.6 


photometric 


Y 




VLT ISAAC 


1980 


2004-11-25 


1.1 


0.7 


clear 


Y 


RX J1856.6-3754 


ISAAC 


1380 


2003-05-22 


1.1 


0.7 


clear 


Y 




ISAAC 


1380 


2003-05-22 


1.1 


0.7 


clear 


Y 


PSR J1932+ 1059 


VLT ISAAC 


1380 


2003-06-05 


1.3 


0.6 


clear 


Y 




VLT ISAAC 


1380 


2003-06-04 


1.3 


0.8 


clear 


Y-mod 




VLT ISAAC 


1380 


2003-06-17 


1.7 


0.7 


clear 


Y 


PSR J2124-3358 


VLT ISAAC 


1380 


2003-06-05 


1.0 


0.5 


clear 


Y 




VLT ISAAC 


1380 


2003-05-22 


1.0 


0.7 


clear 


Y-mod 


RBS 1774 


VLT ISAAC 


1200 


2003-06-06 


1.2 


0.6 


clear 


Y 




VLT ISAAC 


1380 


2003-06-06 


1.2 


0.5 


clear 


Y 


Second Epoch Observations 


Geminga 


VLT ISAAC 


2040 


2007-02-08 


1.5 


0.8 


clear 


Y-mod 




VLT ISAAC 


2040 


2007-02-08 


1.7 


0.7 


clear 


Y-mod 




VLT ISAAC 


2040 


2007-01-26 


1.7 


0.5 


photometric 


Y 




VLT ISAAC 


2040 


2007-01-27 


1.4 


0.7 


photometric 


Y 


RXJ0720.4-3125 


VLT ISAAC 


1980 


2006-12-27 


1.3 


0.8 


clear 


Y 




VLT ISAAC 


3420 


2006-12-27 


1.1 


0.9 


clear 


Y-mod 




VLT ISAAC 


2040 


2006-12-25 


1.1 


0.8 


clear 


Y 




VLT ISAAC 


2040 


2006-12-25 


1.1 


0.9 


clear 


Y 


RXJ1856.6-3754 


NTT 5 OF/ 


24240 


2006-05-09 


1.3 


0.8 


clear 


Y-mod 




NTT S OF I 


15960 


2006-05-10 


1.4 


0.7 


cirrus 


Y-mod 




NTT S OF I 


9120 


2006-05-11 


1.5 


0.8 


cloudy 


Y-mod 


PSR J1932+ 1059 


NTT 5 OF/ 


7200 


2006-05-10 


1.4 


0.7 


cirrus 


Y-mod 




NTT S OF I 


14400 


2006-05-11 


1.6 


0.8 


cloudy 


Y-mod 




VLT NACO 


1360 


2007-05-02 


1.3 


0.6 


clear 


Y-mod 



Appendix C: PSR J2124-3358 

PSR J2124-3358 is a milliseco nd radio p ulsar. The up-to-date 
ATNF pulsar catalogue ( Manchester et al.l2005,) indicates an age 
of 5.86 • lO' years which is considerably older tha n we thought 
when selecting this source some years ago from iTavlor et alj 
d 1 993b in its latest onhne version (no age for PSR J2I24-3358 
noted). At this age and a distance of 270 pc, cooling substellar 
companio ns would have magnit udes fainter than H - 24.9 mag 
following [Burrows et all (Il997l) . Our ISAAC //-band observa- 
tions are not deep enough to detect such objects. We derived 
the NIR limit at the position of the NS which is ///„„ + Icr - 
21.22 + 0.16 mag, the area around the NS is shown in Fig. lC.II 
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Fig. C.l. Detection limit in the //-band for PSR J2124-3358 

The radius for the circle around the position of the NS has a 
radius of 1.5". The box size is roughly 16" x 13". 



